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i-S3 ' Abstract 

Nuclear parton distributions /a(x, Q 2 ) are studied within a framework of the 
DGLAP evolution. Measurements of FA jF® in deep inelastic IA collisions, and 
Drell-Yan dilepton cross sections measured in pA collisions are used as con- 
straints. Also conservation of momentum and baryon number is required. It 
is shown that the calculated Q 2 evolution of F 2 Sn /F 2 c agrees very well with the 
recent NMC data, and that the ratios Rf = /a/ f are only moderately sensitive 
to the choice of a specific modern set of free parton distributions. For general 
use, we offer a numerical parametrization of Rf(x,Q 2 ) for all parton flavours / 
in A > 2, and at 10~ 6 < x < 1 and 2.25 GeV 2 < Q 2 < 10 4 GeV 2 . 



1 Talk given by K.J. Eskola in Quark Matter '99, 12 May, 1999, Torino, Italy. For the transparencies, see 



http://www.qm99.to.infn.it/program/qmprogram.htm] 



1 The Framework 

The motivation to study nuclear parton distributions are the hard probes of strongly inter- 
acting matter [|IJ. Due to a large momentum (or mass) scale Q ^> Aqcd involved, these 
processes take place during the very first fractions of fm/c in the high energy nuclear col- 
lisions, acting as probes of the forming quark gluon plasma. Due to the large scale, hard 
processes are computable within QCD perturbation theory. To a first approximation the 
cross sections of hard processes in nuclear collisions can be factorized as in hadronic colli- 
sions, into parton densities and a hard parton-parton cross section. The parton distributions 
in bound and free protons, however, are different: /U(a;, Q 2 ) ^ f(x,Q 2 ). Typically in hard 
collisions x ~ Q/y/s, so quite different regions in x become relevant when moving from the 
present SPS-energies y/s/A ~ 20 GeV up to the future LHC energies y/s/A ~ 5.5 TeV. 
Thus there is a need for analyses of nuclear parton distributions which consistently cover 
sufficiently wide ranges in x and Q 2 . In this talk, I will present the main results from such 
an analysis ||. ||] . 

In deeply inelastic I A scatterings (DIS), ratios of measured differential cross sections, 
\ dxdQ^ / 1 dxdQ^ ' re fl ec t the corresponding ratios of the nuclear structure function F^ and 
that of deuterium F 2 D . The ratio FA/F® is observed to deviate clearly from unity 0]. Since 
F 2 (x,Q 2 ) = J2 q e l[ X( l( x , Q 2 ) + xq(x, Q 2 )}, parton distributions in bound protons obviously 
differ from those in the free proton. Often the nuclear modifications are referred to as shad- 
owing (x<0.1), anti-shadowing (0.1<x<0.3), EMC effect (0.3 ^x< 0.7) and Fermi motion 
(x — > 1 and beyond). The dependence on the Bjorken x has been observed already in the 
80's [[| but the weaker Q 2 dependence was detected only fairly recently by the NMC ||. 

In Refs. [0, [| our goal has not been to study the actual origin of the modifications but, 
rather - because perturbative QCD (pQCD) does not predict the absolute parton distribu- 
tions - to use the observed effects as input for an analysis in the perturbative region. The 
basic idea in our study is the same as in the global analyses of parton distributions of the free 
proton (like in Ref. ||): we determine the nuclear parton densities at a wide range of x and 
Q > Qo 3> Aqcd through their perturbative QCD (DGLAP 0) evolution by using available 
experimental data and conservation of momentum and baryon number as constraints. 

Information of the nuclear parton distributions can be obtained from I A DIS and Drell- 
Yan (DY) measurements in pA collisions. In these measurements, the accessible values of x 
and Q 2 are strongly correlated, as illustrated in Fig. 1. To perform the DGLAP evolution 
of the parton densities, however, the initial distributions are needed along a fixed scale Qq. 
Therefore, we determine the initial nuclear parton distributions at Qq iteratively, through the 
DGLAP evolution, by using the available data at scales Q 2 > Qq as constraints. Note that 
now the problem is more complicated than in the free proton case because of the additional 
variable A. 

We take the parton distributions of the free proton as accurately known. We choose 
Qq = 2.25 GeV 2 which is the c-quark mass threshold in the set GRVLO ]T2| we are using as 
the basis. We first parametrize the ratio Rp 2 (x, Qq) for (isoscalar) A and x. The potential 
but small nuclear effects in deuterium, and the small tails at x > 1 are neglected here. At the 
initial scale Qq (but only at Qq), we assume that the nuclear sea quarks and antiquarks are 
modified approximately with the same profile, and similarly for the valence quarks. The ratio 



Figure 1: The correlation of x and (Q 2 ) 
in the measurements of DIS || [TT| in I A 
and DY (x = Xq) [JTTJ| in pA. Our choice 



for the initial scale Q% is also indicated. 




FA I F® can then be simply written as a linear combination of a sea quark ratio Rg = Sa/S 
and a valence quark ratio R v = Va/V, where S (V) are the total sea (valence) distributions. 
The ratios Rg and Ry at Ql in turn are constrained by the DIS data ||, |9|, [I(J and DY 
data in pA collisions |]TIJ (see Fig. 1) at higher scales. Ry is also further constrained by 
baryon number conservation. Momentum conservation gives an overall constraint for the 
gluon distributions at Qq. In lack of any direct constraints for the nuclear gluons from the 



Rf 



at very small values of x. The value of x where 
13]. For more details, please see Ref. 



data, we assume that initially 9 a/ 9 

Rq(x, Qq) = 1 is estimated on the basis of Ref. 

0. Once the initial distributions for all parton flavours have been determined like this, the 

DGLAP evolution to larger scales can be performed, and comparison with the data can be 

made (see Fig.l). The initial ratios Rg(x), R v (x), Rq{x) at Ql are then iterated until a 

"best" initial condition is found. 

Fig. 2a shows the scale evolution of the nuclear effects in parton distributions for an 
isoscalar nucleus v4=208. The ratios 9 a/ 9-, Sa/S, Va/V and F^/F® are shown as functions 
of x at fixed values of Q 2 = 2.25 GeV 2 (solid lines), 5.39 GeV 2 (dotted), 14.7 GeV 2 
(dashed), 39.9 GeV 2 (dotted-dashed), 108 GeV 2 (double-dashed), equidistant in logQ 2 , and 
10000 GeV 2 (dashed). For R v only the first and last ones are shown. 

In Fig. 2b we plot the calculated lowest order, leading twist pQCD evolution for the ratio 
F 2 Sn /F 2 c at different fixed values of x, and compare the results directly with the corresponding 
data of NMC ||. The agreement is very good. Note that in order to reduce the potential 
gluon fusion |TJ| effects in the evolution as much as possible, we have chosen the initial scale 
above 1 GeV but below the m c -threshold in order to make the treatment of the initial state 
as simple as possible. The logQ 2 slopes of F 2 , and also of the ratio F^/F®, can be used to 
constrain the gluon distributions fl3| . Unfortunately, however, the values of x of the NMC 
data || are not quite small enough to get a firm handle on the initial nuclear shadowing of 
gluons. 




(b) 




Figure 2: (a) The scale evolution of the ratios xgA/xg, xSa/xS, xVa/xV and F^/F^ for an 
isoscalar nucleus A = 208. (b) The calculated Q 2 dependence of F^/F^ compared with the NMC 
data [Hi. 



2 The EKS98-parametrization 

We have also repeated the analysis by using the CTEQ4L parton distributions || as the basis 
0. Note that in CTEQ4L the sea is more flavour-asymmetric than in GRVLO, and that 
there are less gluons in CTEQ4L at very small values of x. The nuclear ratios Rf = /a// 
for each flavour of partons were, however, found to deviate at most a few per cents relative 
to those computed with GRVLO. We therefore conclude that to a good first approximation 
/a(x, Q 2 )set = Rf (x, Q 2 )f(x, Q 2 )set, where "set" refers to any modern lowest order set of 
parton distributions for the free proton, and where Rf does not depend on the set. 

For practical applications of computing hard cross sections in high energy nuclear col- 
lisions, we have also prepared a parametrization of the nuclear ratios Rf(x,Q 2 ) for each 



parton flavour / in any nucleus A > 2, at 1(T 6 < x < 1 and 2.25 GeV 2 < Q 2 < 10 4 GeV 2 . 
The parametrization is intended for general use, and it is available from us via email, or 
from |http://fpaxpl.usc.es/phenom/| or from [http://www.urhic.phys.jyu.ri/ . 



Finally, we note that our analysis can be improved in obvious ways: more quantitative 
error analysis must be implemented, next-to-leading order DGLAP evolution must be done, 



the initial gluon distributions and gluon fusion corrections |14] must be studied in more 
detail. 
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